Photonic Bloch oscillations of correlated particles 
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A photonic realization of Bloch oscillations (BOs) of two correlated electrons that move on a one-dimensional 
periodic lattice, based on spatial light transport in a square waveguide array with a defect line, is theoretically 
proposed. The signature of correlated BOs, such as frequency doubling of the oscillation frequency induced by 
particle interaction, can be simply visualized by monitoring the spatial path followed by an optical beam that 
excites the array near the defect line. © 2011 Optical Society of America 
OCIS codes: 230.7370, 350.7420, 000.1600 



Bloch oscillations (BOs), i.e., the oscillatory motion of 
electrons in a periodic periodical induced by a dc field, 
are one of the most striking predictions of the semiclas- 
sical theory of electronic transport. BOs manifest the 
wave properties of the electrons, and therefore appear 
in other systems of waves in tilted periodic potentials. 
BOs have been observed for electrons in semiconductor 
super lattices, matter waves in optical lattices, and acous- 
tic or light waves in periodic media. In optics, analogues 
of BOs occur in different structures, including dielectric 
waveguide arrays [1-6], optical superlattices [7,8] and 
metal-dielectric structures [9-11]. Quantum signatures 
of nonclassical light undergoing BOs have been investi- 
gated as well [12-14]. One of the main limitations of pho- 
tonic BOs so far realized is to mimic the motion of single 
particles solely. For many particles, the onset of BOs is 
greatly affected by particle interactions [15-19], and in- 
teresting novel phenomena are predicted for BOs of few 
interacting particles [16,17,19], such as the frequency 
doubling of BOs of two correlated electrons [16, 17]. Ex- 
periments aimed to observe BOs of few correlated elec- 
trons or bosons are rather difficult to be performed in 
semiconductor superlattices or cold atoms, where the 
many particle regime is generally of easier access. In this 
Letter a photonic realization of BOs for two correlated 
electrons is proposed, which is based on light transport 
in a two-dimensional square waveguide array with a de- 
fect line. 

The motion of two interacting electrons moving on 
a one-dimensional tight-binding lattice subjected to an 
external dc force F is described by the Hubbard Hamil- 
tonian (see, for instance, [16,17]) 

H = ~K ( a n+l,s<s +4,s a n+l,*) M 
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where a n?s and a) n s are the annihilation and creation op- 
erators for the electron at site n with spin s =t, I, ft is 
the position operator, n is the hopping amplitude, a is 
the lattice period, and U is the on-site electron-electron 
interaction strength. In order to allow for double occu- 



pancy of the on-site orbital, the two electrons are as- 
sumed to have opposite spins (singlet state). A photonic 
realization of the Hubbard Hamiltonian (1) for two elec- 
trons can be readily obtained after expanding the state 
vector of the system \^(t)) as a superposition of Wan- 
nier states, i.e. \ip(i)) = J2 n ,m c n,m(t)\ns 1 ,ms 2 ), where 
the ket \nsi,ms2) represents a state with one electron 
with spin si at site n and the other electron with spin 
S2 at site m. In the Wannier representation, the time 
evolution of the quantum state ^(t)), governed by the 
Schrodinger equation with h = 1, i.e. idt\ijj{t)) = H\ip(t), 
reads explicitly 
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In their present form, Eqs.(2) can be viewed as the 
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Fig. 1. (Color online) (a) Schematic of a square waveg- 
uide array with a line defect that realizes the Hubbard 
Hamiltonian (1) for two electrons with opposite spins, 
(b) Band structure of the tight-binding lattice for a non- 
vanishing interaction and in the absence of the dc force. 

coupled-mode equations describing light transport in a 
two-dimensional square array of waveguides with a su- 
perimposed transverse index gradient and with a defect 
line (interface) at the diagonal lattice sites n = m, in 
which the hopping amplitude n defines the coupling con- 
stant between adjacent waveguides, the on-site electron 
interaction U defines the propagation constant shift of 
the guided modes for the waveguides along the interface 
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n = ra, and the temporal evolution of the Wannier am- 
plitudes c ny7n is mapped into the spatial evolution of the 
modal amplitudes of light waves, along the array axis, 
trapped in the various waveguides [see Fig. 1(a)]. Such 
a structure could be realized, for instance, by femtosec- 
ond laser writing in fused silica, in which the defect line 
is readily obtained by varying the writing speed of the 
laser beam (see, for instance, [6,20]). The transverse in- 
dex gradient is achieved by circularly-curving the waveg- 
uides in the (Y, z) plane, where the Y direction defines 
the defect line [see Fig. 1(a)] . In the waveguide refer- 
ence frame, propagation of light waves at wavelength A 
is described by the optical Schrodinger equation for the 
electric field envelope <fi(x,y,s) [21] 
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where ft = X/(2tt) is the reduced wavelength of light, s 
is the curvilinear coordinate along the axis of a refer- 
ence waveguide, n s is the substrate refractive index, R 
is the radius of curvature of waveguides, and V(x,y) — 
n s — n(x, y) is the optical potential which is determined 
by the refractive index change n s — n(x, y) inscribed in 
the substrate. Indicating by h(x,y) the normalized in- 
dex profile of each waveguide core and by An, Ani the 
index changes realized in the lattice and along the de- 
fect line, one can write V(x, y) = — Y2 n ml^ n + (^ n i — 
An)S nyrn ]h(x — na, y — ma), where a is the lattice period. 
Coupled- mode equations (2) are obtained from Eq.(3) 
in the tight-binding approximation and neglecting cross 
coupling, provided that time t is replaced by the spatial 
curvilinear distance s and a forcing F = 2nn s / '(\/2R\) 
is assumed. 

To understand the appearance of frequency doubling 
in the BOs of two interacting electrons, it is worth con- 
sidering the energy spectrum of the two-electron Hamil- 
tonian in the absence of the external force, i.e. for F = 
[16,17]. For U = 0, there exists a single tight-binding 
band of energy —4ft < E < 4ft, which is precisely the 
energy band of a single electron on a two-dimensional 
square lattice. The corresponding eigenfunctions are 
Bloch states which are fully delocalized in the crystal. 
For a non- vanishing interaction U, the energy spectrum 
is composed by two bands [16]. The first one, which cov- 
ers the range —4ft < E < 4ft [band 1 in Fig. 1(b)], cor- 
responds to delocalized Bloch states which are scattered 
off by the defect line n = m. The second band [band 2 in 
Fig. 1(b)], which covers the range U < E < \J\J 2 + 16ft 2 , 
is related to the appearance of deject modes which are 
localized at the n = m interface [i.e. around the Y 
axis of Fig. 1(a)]. These two bands starts to separate 
when the interaction U increases above 4ft [see Fig. 1(b)]. 
The initial distribution c nj7n (t = 0) of the two elec- 
trons in the lattice is assumed to be Gaussian-shaped, 
i.e c njrn (t = 0) = Zexp[-(n-n ) 2 /w 2 - (m-m ) 2 /w 2 ], 
were Z is the normalization factor, no and mo are the 
mean positions of the two electrons in the lattice, and 
w measures the localization length (in units of the lat- 



tice period a) of the electronic wave function. Such an 
initial condition is readily realized in the optical system 
by initial excitation of the array, at the 5 = input 
plane, with a Gaussian beam of spot size wa with nor- 
mal incidence at the lattice sites (no, mo), i.e. </>(#, 0) = 
exp[— (x — noa) 2 / (wa) 2 — (y — mod) 2 / (wa) 2 ]. Note that, 
if no = mo , the wave packet describes two electrons that 
initially occupy the same site and thus strongly inter- 
act. Conversely, if | no — mo | is much larger than w, the 
initial wave packet describes two electrons that are ini- 
tially separated each other. When the external force F 
is applied, in the latter case the two particles basically 
undergo independent BOs with a characteristic period 
given by Tg = 2ir/(Fa), provided that the amplitude 
of BOs is smaller than the electronic separation. In the 
optical lattice realization of Fig. 1(a), such a result can 
be simply explained by observing that in this case the 
two-dimensional wave packet motion basically remains 
confined in a homogeneous region of the lattice and does 
not touches the defect line n = m, thus realizing a two- 
dimensional BOs motion. Conversely, if the two electrons 
initially occupy the same site, i.e. for no = mo, BOs 
develop a frequency doubled component which is more 
pronounced at intermediate couplings. As discussed in 
Ref. [17], such frequency doubling is associated with the 
excitation of bounded states of two electrons in a sin- 
glet configuration (i.e. the localized modes at the defect 
line m = n). The relative contributions of fundamen- 
tal and frequency double components of the BOs basi- 
cally depends on the relative excitation of bounded and 
unbounded states from the initial wave packet. In par- 
ticular, the frequency-doubled component is more pro- 
nounced for U ~ 4ft, i.e. at the interaction strength 
at which the two lattice bands starts to separate each 
other, whereas it vanishes in the low and strong coupling 
limits [17]. The onset of independent BOs for spatially- 
separated electrons, and of correlated BOs with fre- 
quency doubling of the oscillation frequency for closely- 
spaced interacting electrons, is shown in Figs. 2 and 3, 
respectively. The figures show results obtained by nu- 
merical simulations of the paraxial wave equation (3) 
in an array of length 5 cm for parameter values A = 980 
nm, n s = 1.522, An = 0.01, Am = 0.01035, a = 8.6 /xm, 
R = (30/ v^) 

cm and for a Gaussian profile h(x, y) of the 
waveguide core of radius 3 fim. The values of coupling 
constant ft, propagation constant shift U and gradient 
parameter Fa entering in the coupled-mode equations 
(2) are estimated to be ft = 3.977 cm -1 , U ~ 4ft, and 
Fa = 2.7973 cm -1 , respectively. The spatial period of 
BOs is Tb = 2it/(Fa) ~ 2.25 cm. Parameter values of 
the input Gaussian beam are no = 10, mo = 0,iu = 1.3 
in Fig. 2, and no = 0, mo = 0,w = 1.3 in Fig. 3. The dot- 
ted curves in Figs. 2(a) and 3(a) show the numerically- 
computed evolution of the beam centroid (n(s)), defined 
by the relation 
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clearly showing the doubling of the oscillation frequency 
when the lattice is initially excited at the defect line 
[compare Figs. 2(a) and 3(a)]. The curves are well approx- 
imated by the corresponding ones for the wave packet 
centroid (n(t)) = ^ n>m (n - n )\c n , m (t)\ 2 obtained by 
numerical simulations of the coupled mode equations 
(2), which are depicted by the solid curves in the fig- 
ures. Snapshots of \<j)(x, y, s)\ at a few propagation dis- 
tances s are also shown in Figs. 2(b) and 3(b). Note that, 
while for two electrons initially far apart each other the 
wave packet undergoes two-dimensional BOs in a homo- 
geneous region of the array, without touching the defect 
line Y [see Fig. 2(b)], for the two interacting electrons the 
wave packet remains strongly localized near the defect 
line, exciting the defect modes at the interface [Fig. 3(b)]. 

In conclusion, a photonic realization of BOs for two 
correlated electrons has been proposed, which is based 
on light transport in a square waveguide array with a 
defect line. Such an optical setting should provide an 
experimentally accessible laboratory tool for the obser- 
vation of interaction-induced frequency doubling of BOs. 
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Fig. 2. (Color online) BOs for two initially-separated 
(no n- interacting) electrons, (a) Evolution of the beam 
centroid versus the normalized propagation length s/Tb, 
and (b) snapshots of \<fi(x,y,s)\ at a few propagation 
distances s. 
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Fig. 3. (Color online) Same as Fig. 2, but for interacting 
electrons. 



4. H. Trompeter, W. Krolikowski, D.N. Neshev, A.S. Desy- 
atnikov, A. A. Sukhorukov, Y.S. Kivshar, T. Pertsch, U. 
Peschel, and F. Lederer, Phys. Rev. Lett. 96, 053903 
(2006). 

5. N. Chiodo, G. Delia Valle, R. Osellame, S. Longhi, G. 
Cerullo, R. Ramponi, P. Laporta, and U. Morgner, Opt. 
Lett. 31, 1651 (2006). 

6. F. Dreisow, A. Szameit, M. Heinrich, T. Pertsch, S. 
Nolte, A. Tiinnermann, and S. Longhi, Phys. Rev. Lett. 
102, 076802 (2009). 

7. R. Sapienza, P. Costantino, D. Wiersma, M. Ghulinyan, 
C.J. Oton, and L. Pavesi, Phys. Rev. Lett. 91, 263902 
(2003). 

8. V. Agarwal, J. A. del Ro, G. Malpuech, M. Zamfirescu, 
A. Kavokin, D. Coquillat, D. Scalbert, M. Vladimirova, 
and B. Gil, Phys. Rev. Lett. 92, 097401 (2004). 

9. W. Lin, Z. Xiang, G. P. Wang, and C. T. Chan, Appl. 
Phys. Lett. 91, 243113 (2007). 

10. A.R. Davoyan, I.V. Shadrivov, A. A. Sukhorukov, and 
Y.S. Kivshar, Appl. Phys. Lett. 94, 161105 (2009). 

11. R.-C. Shm, Y.-C. Lan, and C.-M. Chen, Opt. Lett. 35, 
4012 (2010). 

12. S. Longhi, Phys. Rev. Lett. 101, 193902 (2008). 

13. A. Rai and G.S. Agarwal, Phys. Rev. A 79, 053849 
(2009). 

14. Y. Bromberg, Y. Lahini, and Y Silberberg, Phys. Rev. 
Lett. 105, 263604 (2010). 

15. A. Buchleitner and A.R. Kolovsky, Phys. Rev. Lett. 91, 
253002 (2003). 

16. F. Claro, J. F. Weisz, and S. Curilef, Phys. Rev. B 67, 
193101 (2003). 

17. W. S. Dias, E. M. Nascimento, M. L. Lyra, and F. A. B. 
F. de Moura, Phys. Rev. B 76, 155124 (2007). 

18. J. K. Freericks, Phys. Rev. B 77, 075109 (2008). 

19. R. Khomeriki, D.O. Krimer, M. Haque, and S. Flach, 
Phys. Rev. A 81, 065601 (2010). 

20. A. Szameit and S. Nolte, J. Phys. B 43, 163001 (2010). 

21. S. Longhi, Laser and Photon. Rev. 3, 243 (2009). 



3 



References with full titles 



1. T. Pertsch, P. Dannberg, W. Elflein, A. Brauer, and F. 
Lederer, "Optical Bloch Oscillations in Temperature Tuned 
Waveguide Arrays", Phys. Rev. Lett. 83, 4752 (1999). 

2. R. Morandotti, U. Peschel, J. Aitchinson, H. Eisenberg, 
and Y. Silberberg, "Experimental observation of linear and non- 
linear optical bloch oscillations", Phys. Rev. Lett. 83, 4756 (1999). 

3. H. Trompeter, T. Pertsch, F. Lederer, D. Michaelis, U. 
Streppel, A. Brauer, and U. Peschel, "Visual Observation of Zener 
Tunneling," Phys. Rev. Lett. 96, 023901 (2006). 

4. H. Trompeter, W. Krolikowski, D.N. Neshev, A.S. Desyat- 
nikov, A. A. Sukhorukov, Y.S. Kivshar, T. Pertsch, U. Peschel, 
and F. Lederer, "Bloch Oscillations and Zener Tunneling in 
Two-Dimensional Photonic Lattices", Phys. Rev. Lett. 96, 053903 
(2006). 

5. N. Chiodo, G. Delia Valle, R. Osellame, S. Longhi, G. 
Cerullo, R. Ramponi, P. Laporta, and U. Morgner, "Imaging of 
Bloch oscillations in erbium-doped curved waveguide arrays", 
Opt. Lett. 31, 1651 (2006). 

6. F. Dreisow, A. Szameit, M. Heinrich, T. Pertsch, S. Nolte, A. 
Tiinnermann, and S. Longhi, "Bloch-Zener Oscillations in Binary 
Superlattices" , Phys. Rev. Lett. 102, 076802 (2009). 

7. R. Sapienza, P. Costantino, D. Wiersma, M. Ghulinyan, 
C.J. Oton, and L. Pavesi, "Optical Analogue of Electronic Bloch 
Oscillations", Phys. Rev. Lett. 91, 263902 (2003). 

8. V. Agarwal, J. A. del Ro, G. Malpuech, M. Zamfirescu, 

A. Kavokin, D. Coquillat, D. Scalbert, M. Vladimirova, and 

B. Gil, "Photon Bloch Oscillations in Porous Silicon Optical 
Superlattices", Phys. Rev. Lett. 92, 097401 (2004). 

9. W. Lin, Z. Xiang, G. P. Wang, and C. T. Chan, "Spatial 
Bloch oscillations of plasmons in nanoscale metal waveguide 
arrays ", Appl. Phys. Lett. 91, 243113 (2007). 

10. A.R. Davoyan, I.V. Shadrivov, A. A. Sukhorukov, and Y.S. 
Kivshar, "Plasmonic Bloch oscillations in chirped metal-dielectric 
structures", Appl. Phys. Lett. 94, 161105 (2009). 

11. R.-C. Shiu, Y.-C. Lan, and C.-M. Chen, "Plasmonic Bloch 
oscillations in cylindrical metaldielectric waveguide arrays " , Opt. 
Lett. 35, 4012 (2010). 

12. S. Longhi, "Optical Bloch Oscillations and Zener Tunnel- 
ing with Nonclassical Light", Phys. Rev. Lett. 101, 193902 (2008). 

13. A. Rai and G.S. Agarwal, "Possibility of coherent phenom- 
ena such as Bloch oscillations with single photons via W states" , 
Phys. Rev. A 79, 053849 (2009). 

14. Y. Bromberg, Y. Lahini, and Y Silberberg, "Bloch Oscilla- 
tions of Path-Entangled Photons", Phys. Rev. Lett. 105, 263604 
(2010). 

15. A. Buchleitner and A.R. Kolovsky, "Interaction-induced 
decoherence of atomic Bloch oscillations", Phys. Rev. Lett. 91, 
253002 (2003). 

16. F. Claro, J. F. Weisz, and S. Curilef, "Interaction-induced 
oscillations in correlated electron transport", Phys. Rev. B 67, 
193101 (2003). 

17. W. S. Dias, E. M. Nascimento, M. L. Lyra, and F. A. B. F. de 
Moura, "Frequency doubling of Bloch oscillations for interacting 
electrons in a static electric field", Phys. Rev. B 76, 155124 (2007). 



18. J. K. Freericks, "Quenching Bloch oscillations in a strongly 
correlated material: Nonequilibrium dynamical mean-field theory 
", Phys. Rev. B 77, 075109 (2008). 

19. R. Khomeriki, D.O. Krimer, M. Haque, and S. Flach, 
"Interaction-induced fractional Bloch and tunneling oscillations", 
Phys. Rev. A 81, 065601 (2010). 

20. A. Szameit and S. Nolte, "Discrete optics in femtosecond- 
laser- written photonic structures," J. Phys. B 43, 163001 (2010). 

21. S. Longhi, "Quantum-optical analogies using photonic struc- 
tures", Laser and Photon. Rev. 3, 243 (2009). 



